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We develop an integral form for the bispectrum in general single-field inflation whose domain of 
validity includes models of inflation where the background evolution is not constrained to be slowly 
varying everywhere. Our integral form preserves the squeezed-limit consistency relation, allows for 
fast evaluation of the bispectrum for all triangle configurations expediting the efficient comparison 
of slow-roll violating models with data, and provides complete and compact slow-roll expressions 
correct to first order in slow-roll parameters. Motivated by the recent Planck results, we consider as 
an example a sharp step in the warped-brane tension of DBI inflation and provide analytic solutions 
for the peak of the resulting bispectrum. For the step in the warp that reproduces the oscillations in 
the power spectrum favored by the Planck data, the corresponding equilateral bispectrum is both 
extremely large and highly scale dependent. The bispectrum serves as a means of distinguishing 
such a model from alternative scenarios that generate otherwise indistinguishable power spectra, 
such as a step in the potential in canonical single-field inflation. 



I. INTRODUCTION 

In this paper we develop a technique for calculating 
the bispectrum of primordial fluctuations in general the- 
ories of single-field inflation beyond the usual slow-roll 
assumption where all of the slow-roll parameters are con- 
sidered to be small and constant. This technique is based 
on the generalized slow- roll (GSR) approach [TH5] which 
has been successfully applied to the bispectrum in canon- 
ical single field inflation [THS] . Here we extend this treat- 
ment to also consider terms in the effective field theory 
of inflation associated with time variation of the sound 
speed or equivalently an inflaton Lagrangian that is a 
general function of the field and its kinetic term. Non- 
Gaussianity in such models can be doubly enhanced due 
to a low and variable sound speed. 

This type of enhanced non-Gaussianity can be realized 
in several different ways. In Dirac-Born-Infeld (DBI) in- 
flation [TU1 [IT] , variations in the sound speed are associ- 
ated with features in the warped brane tension [T21 [JJ5] . 
For example a step in the tension produces a step in the 
sound speed. Step features in the sound speed can be 
considered in a more general context such as the effec- 
tive field theory of inflation [IT] . 

Furthermore, transient variation in the sound speed 
has been shown to arise naturally in effective field the- 
ories of inflation from scenarios within multi-field infla- 
tion where heavy fields are integrated out [TH [H] (see 
also [H])- I n these scenarios, turning trajectories in field 
space result in variations in the speed of sound of the 
fluctuations. 

In both of these cases, if the sound speed variations 
are sharp enough they generate characteristic oscillations 
in the curvature power spectrum [6l [18] and bispectrum 
[TBI [TTj] . They represent specific cases of the general phe- 
nomenon that slow-roll violations in single-field inflation 
lead to features in the spectra of curvature fluctuations 
[20H23] . Starobinsky first noted that a sharp change in 
the slope of the inflaton potential lead to oscillatory fea- 



tures in the power spectrum [24. and it has been recently 
shown that these models have large scale-dependent bis- 
pectra [2"5H2"T] . Violation of slow-roll via a rapidly vary- 
ing sinusoidal component of the potential leading to res- 
onance effects in the correlation functions was first noted 
by [221 [23] before it was found to arise naturally in axion- 
monodromy inflation [2"5H3"0] (see also [3DI32j)- Further 
work on slow-roll violating models includes [33H39] . 

On the observational side, oscillations in the curvature 
power spectrum due to transient violations of slow-roll 
were invoked to explain broad glitches in the WMAP cos- 
mic microwave background (CMB) angular power spec- 
trum (4TM2"] . Intriguingly there is also a slight preference 
in the WMAP data for high frequency oscillations in the 
power spectrum near the first peak [8]. In the Planck 
data, this preference persists out to higher multipoles but 
at a somewhat suppressed amplitude and at a harmonic 
of approximately twice the period A£ ss 10 [43] . The 
high frequency oscillations in the power spectrum take 
the same form regardless of whether they came from a 
sharp step in the potential or sound speed and hence it 
does not distinguish between these possibilities [3]. 

Different explanations for features in the power spec- 
trum should be distinguishable in the angular bispectrum 
of the CMB. The techniques developed here allow effi- 
cient computation of the curvature bispectrum in all of 
these cases, from weak violations of the slow-roll approx- 
imation to nearly order unity violations. CMB angular 
bispectrum constraints on oscillating curvature bispec- 
trum shapes were first considered by Ref. [44] using a 
modefunction expansion method [T3] on WMAP data. 
However, the constraints presented there correspond to 
much lower frequency oscillations, for technical reasons. 
A similar analysis was performed on the Planck data [46] , 
where again only periods in the angular bispectrum of 
A£ > 140 were considered and no significant evidence 
for features was found. Neither of these analyses have 
explored the high frequency region of parameter space 
where a large non-Gaussian counterpart of the best fit 



power spectrum feature would be expected. The compu- 
tation of the curvature bispectrum template for match- 
ing power spectrum signatures represents the first step 
for the analysis of high frequency features in the CMB 
bispectrum. 

This paper is organized as follows. In pi} we review 
the in-in formalism and cast the cubic order action in a 
form amenable to bispectrum evaluation. In |III| we em- 
ploy the generalized slow-roll formalism to derive inte- 
gral expressions for the bispectrum that arises in general 
single-field inflation. In §IV| we consider the specific ex- 
ample of warp features in DBI inflation to illustrate our 
technique before concluding in f]V| In Appendix [A] we 
complete the bispectrum description with operators that 
do not appear in DBI inflation as well as treat the re- 
maining slow-roll suppressed contributions. In Appendix 
[E] we use our formalism to derive a compact, complete 
expression for the slow-roll bispectrum to first order in 
slow-roll parameters. In Appendix [Cj we give computa- 
tional details for the DBI example. We take units where 
Mpi = 1/V8ttG = 1 throughout. 



II. BISPECTRUM 

In this section we consider the bispectrum in a general 
single field model of inflation. These models are charac- 
terized by a nearly time-translation invariant expansion 
history through the Hubble parameter H(t) and its as- 
sociated slow-roll parameters 



1 dlnH 

H dt ' 

1 d In en 



2H dt 



(1) 



combined with the sound speed of inflaton fluctuations 
c s (t) and its associated slow-roll parameter 



1 dlnc s 
H dt 



(2) 



Beyond the slow-roll approximation, the parameters ejy, 
T]h and <j\ are allowed to vary with time as long as in- 
flation itself continues without interruption t/j « 1. We 
establish our bispectrum formalism in terms of these gen- 
eral functions that are specified by the model. 

We begin in §11 A| with a brief review of the 'in-in' 
formalism for the calculation of correlation functions and 



in I II B we the show how the cubic action can be written 



in a form that allows a straightforward calculation of all 
bispectrum configurations for models where the slow-roll 
parameters are allowed to evolve. 



mum 



(0(U) = (tf(t*,t )O(t)U(t*,t )), 



(3) 



where U(t*,to) is the time evolution operator in the in- 
teraction picture, 



^(i*,io) = Texp 



i, 

i I Hj(t)dt 

'to 



(4) 



and Hj is the interaction Hamiltonian. We take the ini- 
tial time to to be in the asymptotic past, to — — oo(l-He), 
where the ie prescription projects out the Bunch-Davies 
state initially. We take the final time to be an arbitrary 
epoch during inflation after all of the relevant modes have 
exited the horizon. 

For the curvature bispectrum Bn, we wish to compute 
the correlator 

(£ kl 7t k2 £ k3 ) = (2^) 3 <5(k 1 +k 2 +k 3 )B K (fc 1 ,fc 2 ,fc 3 ), (5) 

where IZ^ is the Fourier transform of the curvature field 
in comoving gauge 



7t(x,£) 



d 3 k 
(2tt) 3 



e lk ' x ft k (i), 



(6) 



and the hat (TV) denotes an operator to differentiate it 
from the mode function (TV). In the interaction picture, 
the curvature fluctuations 1Z evolve according to their 
quadratic action 



S 2 



dtd 3 xC 2 



dtd 3 x^- 



n 2 



% (dJlf 



yielding the equation of motion for 1Z 



d 
dt 



a 3 e H , 



71 - ae H d 2 K = 0. 



(7) 



(8) 



The tree-level bispectrum is then given by expanding 
Eq. (|| [with 0(U) = £ kl (t»)7t k2 (i*)ft k3 (t*)] to linear 
order 



(^k 1 (t.)^k a (*.)^k 3 (*.)) 



2ft 



t* 



<it<£ kl (U )TZk 2 (U )7t k3 (U )Hi(t)) 



(9) 



which can be evaluated in terms of the 2-point correlation 
function or classical modefunctions via Wick's theorem 
by making use of the unequal time correlator 



A. In-in Formalism 

We work in the "in-in" formalism which expresses the 
Appoint correlation function, or more generally, the ex- 
pectation value of a product of field operators 0(t) as 



<ftk(ti)7Mt 2 )) = (27r) J ^(k + k')7e fe (ti)7^(t 2 ). (10) 

Note that we normalize the modefunctions to have di- 
mensions of the square root of the power spectrum 



(ll k (U)Kv(Q) = (2ir) 3 5 3 (k + k')PT l (k) 



(11) 
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We also use the dimensionless power spectrum 



S = 0. In what follows in SIV we shall use DBI to il- 



A 2 R (k) 



2tt 2 



Pn{k), 



(12) 



which is more convenient for expressing dimensionless 
quantities. 

It therefore suffices to calculate the interaction Hamil- 
tonian to cubic order in the curvature perturbation 
[4"5l 05] Hj ~ — / d 3 x£ 3 where the Lagrangian density 
for the curvature field C = C 2 + £3 + ■ ■ ■■ The task of 
computing the bispectrum therefore begins with examin- 
ing the relevant terms in the third order Lagrangian or 
action. 



B. Cubic Action 

In order to compute the bispectrum, we require the cu- 
bic action. We also seek to express its form in a way that 
the relationship between squeezed bispectrum triangles 
and the power spectrum is manifest. 

For a general single field model of inflation, in co- 
moving gauge with gravitational-wave fluctuations set to 
zero, the cubic action is given by |50j 



5, = 



dt d 3 a 



+ °^ L (e H -3 + 3c 2 ) KK 2 



- 2 



KdiFLdiX 



+ 



ae H 



(e H - 2<7i + 1 - c 2 ) K {dny 



?%L± n 2 n 

c{ dt 



d ( a 3 e H e H - r) H 2 • 



dt 



TZ Z TZ 



where 



o X =— s-ft. 



(13) 



(14) 



Here we have written the term usually associated with a 
field redefinition as a boundary term [7, 51 , and dropped 
terms that are suppressed by additional powers of e#. 
Note that we have also dropped boundary terms that 
decay at late times in standard inflation, and thus do 
not contribute to the bispectrum. 

The S term is associated with an operator 1Z 3 that 
does not contribute to squeezed bispectrum configura- 
tions. Its value depends on the specific model of infla- 
tion, and we have maintained generality here in the spirit 
of effective field theory)]] For example in DBI inflation 



1 We do not consider here effective field theory terms associated 
with the extrinsic curvature which appear mainly in ghost infla- 
tion |52| or Galileon interactions 1531 . 



lustrate our technique and so we defer consideration of 
this term to Appendix | A 1| Furthermore the x term only 
contributes at 0(en) to the reduced bispectrum and van- 
ishes for squeezed configurations. For completeness, we 
consider its effect in Appendix |A 2| In the remainder of 
this section, we therefore drop these two terms as neither 
play a role in establishing the consistency of the squeezed 
bispectrum and power spectrum spectral index. 

W hile the remaining four terms in the action of Eq. 
(13) are otherwise complete, each contains c~ 2 enhanced 
terms. In particular the 7?. 2 7?. terms do, and they con- 
tribute to squeezed bispectrum configurations. On the 
other hand, we know that the theory satisfies the consis- 
tency relation 



12 = B K (k s ,k L ,k L ) 
5 lNL ~ ks™o P K (k s )P n (k L ) 
= d\nk 3 L P n (k L ) ^ l 
din 



(15) 



and so cannot have cj 2 enhanced terms that contribute 
to squeezed triangles. While it is well known that the 
consistency relation is satisfied in slow roll by the cancel- 
lation of terms, beyond slow-roll it is difficult to establish 
in this form. 

Our strategy is to combine the c~ 2 terms that form to- 
tal derivatives of quantities which vanish outside of the 
horizon and hence provide no contribution to the bispec- 
trum. We begin by using the equation of motion (Jsj) and 
dropping a total space derivative to show 



2F nc 2 = F^n 2 n) - F^n 2 n, (16) 

at \ ci I ci 



for an arbitrary function of time F(t). Hence with F — 

(e H -Vh)/c 2 



S 3 = / dtd 3 xTZ 



£H -Tf H - a i/ 2 r 

^ n i-2 



Here % 2 is the quadratic Hamiltonian density 

„2 



n 2 



a 3 e H 



ii 2 + ^ [dnf 



(17) 



(18) 



and the quadratic Lagrangian density was given in 
Eq. (0. 

Next we note that several terms in S3 can be grouped 
into a total derivative generalizing Ref. [M] 

-2 (e H -Vh-y) nL ^ 

+ {e H +a 1 + —)nn 2 , (19) 
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where again F is an arbitrary function of time. 
F = l/cl 



With 



S3 



dtdrx 



dt \ Hci 



K{H 2 + 2£ 2 ) 



C-2 



(20) 



The total derivative is irrelevant here, and we can drop 
itj^j Although this form for S3 efficiently groups the terms 
into cj 2 enhanced an non-enhanced terms, it still ob- 
scures the consistency relation from the latter by not 
manifestly scaling with slow-roll parameters. To expose 
this relation, we can reverse the above operations on the 
'H2+2C2 term by using the identities Eq. ( [l9| with F = 1 
and ( 16 1 with F — 2ch — T)n + ci/2 to obtain 



S, 




where we have again dropped the total derivative term for 
the same reasons as above. Note that for canonical fields 
c s = 1 and this form reproduces the r)n dependence of 
the cubic action used in Ref. [7]. It furthermore restores 
6h terms which typically do not contribute significantly 
to the bispectrum but complete the consistency relation 
in the slow-roll limit. 

In fact the slow-roll consistency relation can be triv- 
ially demonstrated given this form for S3 by making use 
of the technique of Ref. 54J. Notice that, to leading or- 
der in slow roll, the only term that contributes to the 
bispectrum in the squeezed limit is the boundary term, 
the second line of in Eq. (21 1. The first term on the last 



line can be seen to be higher order after making use of 
the identity at Eq. (19 1. Since the boundary term is a 



total derivative, we can easily evaluate its contribution 
to the bispectrum. Making use of the Hcrmiticity of the 
fields, we can rewrite Eq. (J9j) as the commutator 



n 



d 3 qj a 3 e H 



(2e ff 



2 



X (27r) 3 <5 3 (qi+q 2 +q 3 ) + ... 



(22) 



2 Recall that H2 generates time translations for 1Z, and so this 
operator results in terms which involve lZ(t t ) which is exponen- 
tially decaying at late times when all modes of interest are outside 
the horizon (at least in standard inflationary scenarios). Since 
the resulting correlation functions vanish exponentially fast at 
late times, they are irrelevant as long as we only interested in 
correlations late in inflation. 







Operator 
V}1l 



Source Squeezed Tij Eq. 



1 TZ(H 2 + 2£ 2 ) 

2 TZC 2 

3 ii 3 

4 TZ (dTZ) d X 



2e H ~r)H + 



yes 
yes 



(421 



Cs 

aHs 

, C a 

aHs 



(461 



(501 (54i 



no (All 1 (A17) 



(A22) 



TABLE I. GSR bispectrum operators, sources, and triangle 
weights. 



1.) 



00, 10, 20, 30, 40 xsinx 

23,33 x sin x + cos x 

01,11,21,22,24,31,32,41 cosx 



02, 12 

25, 34 

26, 35 



- sin a; 



: sin x — cos x 
W{x/2) 



TABLE II. GSR window functions for various operators i and 
a:-we igh t factors with W as the power spectrum window of 
Eq. |38l). 



where '. . .' refers to terms that vanish in the squeezed 
limit. We can then make use of the fact that the canonical 
momenta of this theory is given by, 



dn 2 ah H . 

TTk = -7^- = 2 — s— /c k , 
ci 



and satisfies the canonical commutation relation 

[ft k ,7T k ,] =z(2^) 3 3 (k + k'), 



(23) 



(24) 



to evaluate the commutator in Eq. (22 1. The leading 



order squeezed limit is therefore given by 

B n (ki,k 2 ,k 3 ) w (2e„ -r] H + y) Pn{ki)Pn{k 2 ) 

+ perm. (25) 

where "perm." refers to the two cyclic permutations of 
the hi indices. Since in slow roll 

d In A2, 



n, - 1 = 



R 



d\nk 

= -(4£ H -2riH + <ri), (26) 

and Pn(k 2 )PTz{k3) < Pn{ki)Pn{k 2 ) for k x < k 2 ps 
k 3 , this establishes the slow-roll consistency relation, 
Eq. (pi 



III. GENERALIZED SLOW ROLL 

In this section we construct an efficient integral formu- 
lation of the bispectrum contributions from Eq. ([9]) and 



4 



T% i 



Equilateral 



Squeezed 



Too 


-1 


-1 


Toi 


—6 


k s 


T02 


9 




T w 


7 
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Tn 


17 
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2^ 


Tia 


9 
2 





T20 

T21 

T22 
T23 
T24 
T 25 

T26 

T 2 b + perm. 

T30 
Tai 
T32 
T33 
T 34 

T35 

T AB + perm 



_ 5 
6 

_7 
6 

5 
(i 
1 

4 

_37 
4 

II 

3 

_1S 
2 



-11 

16 
-6 



(ft)' 







1 



III' 







3 
2 

"fcs 2 



+ 6^+3 



Tin 

r 4 i 



TABLE III. GSR triangle weights for equilateral and squeezed 
configurations +0(ks /fez,). Squeez ed c ontrib ution s from T2 m 
and T3 m , m > 3 cancel due to Eq. (571 and (A18l. 



the cubic action of Eq. (21 1 that allows arbitrary time 



variation in the slow- roll parameters €h,i]h,^i- The re- 
maining model-dependent 1Z 3 term and tu suppressed 
terms are considered in Appendix |A 1| and |A 2| respec- 
tively. 



A. Formalism 

To evaluate the bispectrum exactly, we need to solve 
for the background evolution for the slow-roll parameters 
ejj, i]Hi &i aud sound speed c s as well as the modefunc- 
tions IZk for each mode in the triangle configuration. The 
modefunctions themselves are dependent on the slow- 
roll parameters and beyond the slow-roll approximation, 
where the slow-roll parameters are taken to be constant, 
there is no general analytic solution for their behavior. 

The generalized slow-roll approach (GSR) [TJ[51|3] pro- 
vides an iterative approximation to the modefunctions. 



Their equation of motion Eq. ^ can be recast as [5] 

2 



fy 

dx 2 



1 - 



y 



g(lns) 
x i 



where 



2tt 2 x 

x = fcs, and the sound horizon 

ft, 



s(t) 



c s dt 
a 



with t en d defining the end of inflation. Here 

/" - 3/' 



9 



f 



with ' = d/dlns throughout and 

2 €HCs ( aHs 



8ir z 



H 2 



(27) 



(28) 



(29) 



(30) 



(31) 



f- 2 (sec Eq. B8 I 



Note that in the slow-roll limit A^ 
In the GSR approximation, one first defines the solution 
to Eq. (27) with g = and Bunch-Davies initial condi- 
tions 



Vo(x) 



1 



(32) 



and then replaces the RHS of Eq. ( 27 ) with y — > y . The 
solution to first order in g is 



du 



g(lns)y (u)%[ya(u)y (x)], (33) 



where u = ks. With these relations we can define an 
integral approximation to the bispectrum to leading or- 
der in the slow-roll deviations g, e#, rjn and cr%. For £3 
operators that already include slow-roll parameters only 
the zeroth order yo modefunctions are required whereas 
those that involve none, i.e. 1ZC2 for DBI inflation, the 
first-order modefunction correction contribute to first or- 
der in the GSR approximation. 



B. Integral Form 

All contributions from the S3 operators can be cast 
into integral form for the dimensionless bispectrum 



k X k 2 k Z (2tt) 4 A2 



k 2 k 2 k 2 



(34) 



where A s is a constant of order the dimensionless power 
spectrum = fc 3 Pr / '2tt 2 . In the leading order GSR ap- 
proximation these integrals depend only on the perimeter 
of the triangle K = ki + k% + k^ rather than its shape. 
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Thus this description enables a highly efficient compu- 
tation of all bispectrum triangles from a handful of one 
dimensional integrals. 

We will group our integral results according to sources 
Sij (In s) indexed with: (z) the operators they correspond 
to; (j) the scale x = Ks at which the operator sources 
contribute. The integrals are given by 



Iij(K) = SyQnsJWijiKa*) + 



ds 



S' ij (lns)W ij {Ks) 
(35) 

where Wij are fixed window functions that are indepen- 
dent of the source. The triangle shape dependence is 
carried by which are universal functions of (k\, k 2 ,k 3 ) 
such that any bispectrum triangle can be computed as 



Q 



A K (fc 1 )A K (fc 2 )A K (fc 3 ) j ^ 



[T 2B J 26 (2fc 3 )+perm.]}. 



(36) 



In the following sections, we consider the contribution 
of each operator composing S3 in turn. We summarize 
these terms in Tab.[TJ their associated windows in Tab.[TTl 
and limiting cases of the triangle weights in Tab. |III| 

In the leading order GSR approximation, the power 
spectrum itself is given by the integral [5] 



lnA^ = G(lns»)W(fca*) 



ds 



W(ks)G'{lns), (37) 



where the power spectrum window function is 
3sin(2u) 3cos(2u) 3sin(2u) 



and the source function 



2 it 



G = -21n/+-(ln/)'. 



(38) 



(39) 



This power spectrum expression then completes the form 
of the bispectrum in Eq. ( 36 1 . 



0: Wll 



and the triangle weights 



T, 



on 



T01 = - 



i& kik ) 



T02 — 



kik 2 k 3 
hk 2 k 3 



(42) 



are identical to the canonical case. Compared with the 
treatment in Ref. [7], here we include the eu cor rection 
associated with the full TZ 2 TZ operator in Eq. (21). 



As in the canonical case, the appearance of 1/f in 
the source makes the integrals in Eq. (35 1 involve total 



derivatives S[a and hence guarantees that the bispectrum 
remains constant once all 3 fc-modes have exited the hori- 
zon. It was shown in Ref. [7J that this is the consequence 
of first-order modefunction corrections using Eq. ( 33 ) . 

Note that both the T 01 and T 02 terms contribute to 
squeezed triangles k± — ks <C &l = k 2 ss k 3 



lim T 02 



-2T 01 



kr 



(43) 



and hence are involved in establishing the consistency 
relation. 



2. 



1: 1Z{H2 + 2C 2 ) 



The 1Z(H 2 + 2C 2 ) can likewise be calculated with the 
GSR expansion. Again only the zeroth order modcfunc- 
tions are required throughout and constancy of the bis- 
pectrum on superhorizon scales is automatic due to an 
integration by parts which brings the source into the form 



'10 



Sn — S\ 2 



(44) 



with windows 



Wio(x) = a; sin a;, 
Wn(x) = cosx, 

W 12 {x) = , 



(45) 



We start by considering the first two terms in Eq. ( 21 ) 
which involve the 1Z 2 1Z operator. This operator was con- 
sidered in detail in Ref. [7] for canonical (c s = 1) scalar 
fields and the results carry directly over to the general 
case with the source replacement 



"5*00 — "5*01 — S 02 — j (^e H — rj H + 

Namely the windows 

Woo(x) = xsinx, 
W i(x) = cosx, 



W 02 {x) 



SIM 



(40) 



(41) 



and triangle weights 



K 2 



2kik 2 k 3 



'c 3 

i>j i^j 



(46) 



Both Tn and T 22 contribute in the squeezed limit and in 
fact 



kh 

lim Tn = -T12 = 2 — . 

fcs<fei ks 



(47) 
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In the slow-roll limit only the boundary terms in In and 
I12 contribute since source derivatives involve evolution 
in the slow-roll parameters (see Appendix [B] for further 
discussion). Given that V^n(O) = Wi 2 (0) = 1, the sum 
of the two terms vanishes and the operator does not con- 
tribute to squeezed triangles in the slow-roll limit. They 
do contribute beyond the slow-roll limit where the evo- 
lution of the sources enters. 



S. 



2: TZC 2 



The only operator that is enhanced by c~ 2 is 1ZC2- 
This one is also special in that it is not suppressed by 
any slow-roll parameters. To calculate its bispectrum 
consistently to leading order in the GSR approximation, 
we must use the first order modefunction expansion from 
Eq. ( 33 ) [see { A 1 for more details on an analogous calcu- 
lation^ These terms contribute comparably to those that 
describe the evolution of the c7 2 enhanced term itself 



S 2 (lna) 



\aHs 



(48) 



using zeroth order modefunctions uq. 

For the latter effect, it is computationally advanta- 
geous to isolate the S'2 evolution terms by integrating 
by parts so that the new sources are 

5*20 = 5*21 = S' 2 , 

S22 = S 2 , (49) 



with 



T20 — 



[21 



2 X)»>j k i k j 



kik 2 k; 



2K 2 



(50) 



k{ kj 



K ^ 

i>3 



4 

K 2 



1 



kik 2 k 3 12 t-r 1 



i>j 



K 2 X] k i fc i 



and 



W2o(x) = x sin a;, 

W 2 i{x) = W 2 2(x) = cosx. (51) 

By constructing the integrals in this manner, we guaran- 
tee that evolution in S'2 that has compact support in S' 2 ' 
results in rapidly convergent integrals as x — > 00 which 
do not require regulation. This should be compared with 
the unmanipulated S 2 integrals or the exact integration 
over modefunctions, both of which have w indo w or mod- 
efunction weights that diverge as x 3 (see { Al). 



Note that the boundary term of the I 2 2 integral gives 
the well known result that 



1 



- 1 



kf- 



i>j 



1 



2K 2 



k itf 



to zeroth order in slow-roll parameters. In Appendix [B] 
we derive an expression that is correct to first order in 
the slow-roll approximation. 

Unlike the other operators, there are leading order con- 
tributions associated with the deviation of the modefunc- 
tions from their yo de Sitter form. The nested integration 
of the modefunction correction inside the bispectrum in- 
tegral can be unwound by integration by parts, using the 
leading order approximation of S2=const., leaving a new 
source 



9 "5*2 — S 23 — S 24 — S25 — S 2 q, 



(53) 



with 

T23 
T24 

T25 



3 2 J2t>j kikj 
2 ~ K 2 
(K -2h)(K 



2k 2 )(K-2k 3 ) 



Akik 2 k 3 
1 



- 6fcifc 2 4- ^2)^3 



To, 



8kik 2 k 3 (K - 2k 3 )K 2 
x ( fc i - k l) 2 { k l +6fcifc 2 - 
+A(k 1 -k 2 ) 2 (k 1 +k 2 ){k 2 1 
+2(3kf + 2ik\k 2 + 64fc 2 fc| + 23fc x fc| + 3fc^)fc 2 
+16k 1 k 2 (k 1 + k 2 )kl - (7k 2 + 20fcifc 2 + 7k\)k\ 

-4(fci + k 2 )k 3 +perm., 
1 



\2k x k 2 k 3 (K -2k 3 ) 2 K 2 



(h - k 2 ) 2 (k 1 + k 2 f(k 2 + 6fc 1 fc 2 + k 2 ) 



+3(kf - KY{k{ + 6fci/c 2 + k z 2 )k 3 
+2{ki + k 2 )(6kf + 35kfk 2 + 106fc 2 fc 2 



+35fcifc| + Qk^kl 
+2(2fcf + hk\k 2 - 2<6k\k\ 



K)k \ k^ ~\~ 2/^2 



-(fc 1 + fc 2 )(19fc 1 

'2 



^k Y k 2 



I9k 2 )ki 



+(-9k( + Akxk 2 - 9kj)kl 



-&{ k l+ k 2) k \ 



2k 3 



perm. 



(54) 



Here "perm" means the 2 additional cyclic permutations. 
Integration by parts on the modefunction expansion also 
leaves a boundary term that is described by 



1 



L2B 



6hk 2 k 3 {K - 2k 3 ) 



(55) 



2(fci 
(fci- 



k 2 ) 2 (k 1 



+ k 2 )(k 2 1 
2 {k 2 + 3k 1 k 2 - 



- 2>k\k2 
k 2 2 )k 3 



fe)(3fc 2 + 5hk 2 + ikl)k 2 3 
(3fc 2 - 2fcifc 2 + 3fc 2 )/c| 



•2(fci 



k 2 )ki - k 5 



3 ■ 



(52) and 726(2^) in Eq. (36). The windows associated with 
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these terms are 



W23(x) =a;sinx + cosx, 
W 2 ±{x) =cosx, 
sin x 

W 2 $(x) =2 cos a;, 



W 26 (x) =12 



(56) 



4.T 



Note that W 2 &{x) = W(x/2), the window function of the 
GSR power spectrum defined in Eq. (38). 

While the evaluation at in Eq. (35 1 for these 5 mod- 



cfunction terms would formally require integrating gS 2 , 
the sum exactly vanishes for all triangles and Ks* <C 1. 
They may be omitted in practice so long as the same s* 
is taken for each. This cancellation is a consequence of 
the triangle weights obeying 



J=3 



[ 2j 



(T 2 b + perm.) = 0. 



(57) 



Likewise, while individual terms would seem to con- 
tribute to squeezed triangles (see Tab. 
suppressed by kgjkh for any source evo 
that source contributes when k$s <C 1. 



nil, the sum is 



ution so long as 
Note that the 
consistency relation is not expected to hold if the long 
wavelength mode, fcg, is inside the horizon since it can 
no longer be considered as a change in the background 
for the evolution of the short-wavelength, fcj,, modes. 



C. Consistency Relation 

With these integral expressions for the bispectrum and 
power spectrum, we can now re-examine the consistency 
relation for squeezed bispectra, Eq. (15). Amongst the 



integrals only Ioi, I02, hi, III contribute to squeezed 
triangles as k^/ks- Note that none of these terms are 
enhanced by cj 2 relative to the power spectrum. To 
leading order, 



12 



NL 



-2 



/* 



C-H 
f 



Vh 
f 



W v (k L s) + 



W e (k L s) (58) 
W a {k L s) , 



where 



W t {x) = -2cos(2x) + - sin(2x), 



This should be compared with the local slope of the 
power spectrum [9] 



din 



din k 



= 2 



r°° ds 

/ — W'{k L s)G'(\ns) 

Js, 8 

ds (f'\ 



(61) 



W n {k L s), 



where 



WJx) = -2 cos(2a;) + - sin(2a;). 

x 



(62) 



The boundary term obviously matches between 
Eqs. (58) and (61) and establishes the slow-roll consis- 



tency relation. The integral piece contributes when there 
are features that violate the slow-roll approximation. For 
a sharp feature at k^s 3> 1, the parameters with the high- 
est number of derivatives of H and c s dominate and we 
can approximate 



Vh + gi/2 
/ 



(63) 



which matches the t\h term in Eq. ( 58 1 and the o\ term 



in the fc^s > 1 limit assumed. Note that <j\ enters with 
opposite sign relative to r\n between the slow-roll and 
sharp feature expressions. This is a consequence of the 
H 2 + 2C 2 term entering the latter but not the former. 
For k^s ~ 1, the power spectrum source G' no longer ap- 
pears as a sharp function compared to the windows and 
so other terms that impact its shape matter [6] . We illus- 
trate below that the GSR approximation maintains the 
consistency relation even in this region. Finally note that 
once all the terms are considered, including window func- 
tion expansions, the contribution of terms not involved 
in the consistency relation is suppressed by 0(kg /&l) 2 
independently of the sources as expected [53] . 



IV. DBI STEP FEATURE 

In this section, we illustrate the GSR integral construc- 
tion of the bispectrum from |III| in a DBI model with 
sharp features in the sound speed. We review the DBI 
model in §IV A[ test the GSR approximation in §IVB| 
and discuss analytic scaling results in §IV C| Appendix 
[C] gives details on how we set DBI parameters that are 
matched to the Planck data. 



WJx) — 2 cos(2x) sin(2a;), 

x 

W a (x) = cos(2x). 



(59) 



Here we have again evaluated the boundary term by as- 
suming is an epoch during slow-roll 



-(4e H - 2riH + cri)\ a ,. 



(60) 



A. Model 

The DBI action is a specific incarnation of the general 
fc-inflation action 1551 



S = Id xy/—g 



R 



P(X, 



(64) 
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Here R is the Ricci scalar. The scalar field Lagrangian 
is taken to be a general function of the field value <p and 
its kinetic term 



X 



1. 



(65) 



The scalar field behaves as a perfect fluid with pressure 
P, 



2XP X -P, c 2 s =P iX /p, x . 



(66) 



In these models the 2 term in the cubic action Eq. ( 13 ) 
is given by 



s = 4-i-2 A 



where 



H 2 e H 



A = X 2 P^xx 



E' 



X A P 



XXX- 



(67) 



(68) 



In the DBI case, 



P(X, <f>)=[l- y/i - 2X/T(</>)\ T{$) - ^(0), (69) 

where T(<p) gives the warped brane tension, and V(4>) is 
the interaction potential. Note that 3 = and so the 
only terms involved in the bispectrum are those given in 

m 

We illustrate our bispectrum technique on models 
where T(0) has a step feature [13]. Details of this model 
including the background evolution and construction of 
the slow-roll parameters can be found in Ref. [BJ; we re- 
view its basic properties below. The warp factor 



as well as the sound speed at a sound horizon of s — 3.692 
Gpc in the absence of the step (b = 0) as described in 
Appendix [Cj The observables are weakly dependent on 
Vq at fixed s and so we follow Ref. [BJ in choosing 



V Q = 7.10 x 10 



-26 



(74) 



In the presence of a step, the power spectrum and sound 
speed vary from these values but we consider A s as a fixed 
number and label models with c, evaluated at b — > 0. 



The value of the sound horizon in Eq. ( 73 ) is motivated 
by the recent Planck results [33] where a step feature at 
this scale 

s s = s(4> s )= 3.692 Gpc (75) 
with power spectrum amplitude in the DBI model of [6 
_ 1 — c. b 



C 



l + c s y/T^Zb 



0.1. 



(76) 



Unlike the WMAP data [5], the Planck data reveals a 
preference for a finite width with their increased angular 
resolution 143 



Xd 



nd 



87.4 



(77) 



with this set of parameters favored at A2 In C « 12. We 
shall see that by preferring a specific x<i, the Planck data 
favor a specific and large maximum for the bispectrum. 
We therefore illustrate bispectrum results with these val- 
ues for b, 4> s and d below while also testing for robustness 
to parameter variations. 



GSR Tests 



r(0) = -^[i + 6F(^)] 

As 



has a tanh step-like feature 



F{4>) = tanh 



and 4> inflates on a potential 



-1, 



V{4>) = v ( l - -H 2 



(70) 



(71) 



(72) 



rolling from small to large values. We have chosen a 
convention that after the feature, T(4>) goes back to its 
b = value. We assume that inflation ends when 4> = 
0cnd- Note that the step is 2b in amplitude. 

The parameters As, /3, and </> e nd are chosen to fix the 
amplitude and tilt of the power spectrum to the maxi- 
mum likelihood of the Planck data 1561 



A s = 2.69 x 10" 9 , 
- 1 = -0.0381, 



To test the accuracy of the GSR approximation, we 
compare integral approximation at Eq. ( 36 ) to a numer- 



(73) 



ical computation of the full bispectrum using mode-func- 
tions obtained from numerically solving the linear equa- 
tions of motion. 

We begin by testing the squeezed bispectrum limit 
which mainly checks the validity of the method. We 
have constructed the GSR approximation in a manner 
in which the squeezed limit and its relationship to the lo- 
cal slope of the power spectrum is manifest in the i = 0, 1 
terms (sec S III C I. We show in Fig.[l]that the consistency 
relation itself is satisfied in the exact calculation of the 
bispectrum and power spectrum. We choose to fix the 
ratio fcs/fci = 10 . Since kss s <C 1 even for k^ in the 
damping tail, agreement with the consistency relation is 
expected. We know analytically that the i = 2 opera- 
tor contribution is suppressed by {ks/k^) 2 compared to 
true squeezed contributions in this case and so it is com- 
putationally most efficient to drop them outright. Note 
that the full extent of this suppression is achieved from 
cancellation of integral terms which is difficult to repro- 
duce numerically. Nonetheless our numerical integration 




FIG. 1. Test of the consistency relation: n a — 1 computed 
exactly from the power spectrum vs /nl computed exactly 
from the bispectrum for a large amplitude step b = —0.2, 
x d = 15, c s = 0.15, k s /k L = 10~ 4 . 



is sufficient to make residual contributions from numer- 
ical errors negligible here compared with the true ones 
from i — 0,1. 

The accuracy of the GSR approximation for these 
highly squeezed triangles is shown in Fig. [2] In this 
case the approximation has small but notable amplitude 
errors even in the small amplitude case. As shown in 
Ref. [S], these errors in amplitude arise due to slow-roll 
corrections in the value of / between when the mode 
kg left the horizon and when the features at are im- 
printed. For the step feature the latter is fixed at s s and 
the correction is [S] 



R = 1 



n s - 1 . ( k s s s 
— ^ — m 



, k s s s < Xf (78) 



e 2 "^/2 « 2.07 is 



and R — 1 otherwise. Here Xf = e 2 ~ lE /2 
the freezeout epoch for the tilt (see Tab. IV). In Fig. [2] 
we also show that with this correction, the remaining 
error from the feature is in a small out-of-phase compo- 
nent. These too can be corrected with the first order 
techniques of Ref. [9] but note that in this example they 
peak at an unobservable below 0.02. There is an 

additional correction for the slow-roll contributions on ei- 
ther side of the features due to the evolution in / between 
horizon crossing of k$ and k^. Since these are slow-roll 
suppressed, correcting them is never relevant. 

Next we compare the equilateral bispectrum where the 
c~ 2 enhanced i = 2 term contributes. Here we fix the 
basic parameters as described in the previous section ac- 
cording the Planck best fit, but allow the sound speed 
and amplitude b of the feature to vary. 



FIG. 2. Exact vs GSR approximation for squeezed triangles 
for a small amplitude sharp step b — —0.01, x d = 15, c s = 
0.15, ks/kh = 10 -4 . Top panel shows the GSR result with 
a small correction due to ks exit ing the horizon many efolds 
before the feature from Eq. (78 1. Bottom panel shows the 



error in the corrected GSR approximation as well as that of 
the original GSR approximation. After correction, the errors 
are mainly an out of phase component. 



Fig. [3] (top) shows the result for a small amplitude 
feature with b = —0.005 for two different values of 
c s = 0.1,0.5. The GSR approximation captures the am- 
plitude of the features to a few percent for both cases. 
In particular the amplitude increases as (k cq s s ) 2 until it 
reaches a peak at k cq s s ~ 2xd before damping away. As 
we shall see below, varying the sound speed changes the 
relative weight of the different operators and so the two 
different cases demonstrate that the relevant terms have 
been individually calculated correctly. 

In Fig. [3] (bottom) , we test the larger amplitude fea- 
ture preferred by the Planck data. Given that the same 
fractional change in the warp or sound speed causes a 
smaller effect as c s — > 1, the amplitude of the step b re- 
quired to match the data increases with c s . On the other 
hand the accuracy of the GSR approximation depends 
directly on b and so errors increase with b and c s . 

For c s = 0.15, the Planck parameters predict an ex- 
tremely large equilateral bispectrum at peak. We illus- 
trate this value since the maximum violates even a weak 
criteria for the validity of perturbation theory 



1,3 



Ak < 1. 



(79) 



Beyond this point, the curvature field is strongly non- 
Gaussian and loop corrections likely invalidate the cal- 
culation. Note that this sets a firm lower bound on c s 
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FIG. 3. Exact vs. GSR approximation for equilateral triangles. Shown are a small b = —0.005 amplitude warp step (top) and 
the large amplitude step that matches the Planck feature (bottom) for a low c s = 0.15 (left) and high c 3 = 0.5 (right) sound 
speed. For the small amplitude step the GSR approximation is accurate in amplitude at the several percent level. For the large 
amplitude step the amplitude accuracy is ~ 20% or better at peak. Other parameters have been set to reproduce the Planck 
best fit, in particular Xd = 87.4. 



for models that seek to explain the Planck feature with 
steps in the warp. 

For the case of c s = 0.5, the full change across the 
step approaches order unity, specifically \2b\ ~ 0.36 and 
the GSR approximation at the peak of the bispectrum 
holds to ~ 20%. For larger values of c s it is impossi- 
ble to explain the feature without a very large step in 
the warp. Note a step in the potential does not suffer 
this problem and they remain viable explanations even 
as c s — > 1. Likewise for a potential step in canonical sin- 
gle field inflation with the same parameters that fit the 
Planck power spectrum |£//fc^ q | max ~ 766. The ampli- 
tude of the equilateral bispectrum clearly distinguishes 



the two scenarios for c s < 0.5. 

The errors in the GSR approximation at high b are 
mainly due to the i = 2 operator. Note that in obtain- 
ing the modefunction correction source to this term in 
Eq. (53), we have assumed that S2 in Eq. (48) is con- 



stant whereas it is a function of c s . Corrections to this 
approximation from the sound speed step are expected 
to be important as the amplitude of the step increases. 

In fact the errors in the approximation around k cq s s ~ 
1 can be directly attributed to this problem. If X4 3> 1 
subleading terms in the approximation begin to dom- 
inate here and cause order unity errors. Fortunately, 
this occurs only in the region where the bispectrum is 
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FIG. 4. Modified GSR source correction for the highest am- 
plitude b case of Fig. [3] By changing the i = 2 modefunction 
correction source to a total derivative of a slow-roll suppressed 
quantity through Eq. (ISO]) , the errors at ks B ~ 1 decrease. 



too small to be observed due to cosmic variance. We 
can nonetheless correct for this problem by modifying 
the source S2m f° r m > 3 from the S' 2m = gS 2 form 
of Eq. ( 53 1 . The problem is that this form is not a total 



derivative of a slow-roll suppressed source and so for large 
amplitude features it does not integrate back to a slow- 
roll suppressed S2m after the inflaton has transited the 
feature. Analogous effects in the i — operator are fixed 
by self consistently expanding to the next order in the 
GSR approximation [7j. Since this would involve nested 
integrals, we can fix the problem by simply replacing the 
source with a form that is identical in the slow-roll and 
small feature limits but which carries the total derivative 
structure 



S 2m ^^[G(\ns)-G]S 2 



(80) 



To ensure that the source is slow-roll suppressed after 
the feature we set the constant G — G(\ns s ;b = 0). In 
Fig. [4] we show that this change corrects the ks s ~ 1 
problem. On the other hand, this problem appears in a 
non-observable part of the spectrum for this model and 
also disappears if x^ < 25, where one might otherwise 
think ks s ~ 1 effects are important, and so we do not 
consider it further. 



C. Analytic Scaling 

Given the sensitivity of the maximum equilateral bis- 
pectrum amplitude to the choice of parameters and the 
large current observational uncertainties on their values 



FIG. 5. Comparison between exact, GSR and leading order 
analytic solution predictions for the maximum amplitude of 
|5/fce q [. The step parameters are b = —0.005, Xd = 30.2. 
Dashed lines indicate where the individual contributions have 
the opposite sign to total. For c s < 1 the i = 2 term domi- 
nates whereas for c s — > 1 all terms are comparable. 



from the Planck data, it is useful to have analytic scalings 
for the amplitude and form around maximum. 

In the sharp-feature limit, the dominant contributions 
to the bispectrum arise from the terms with the highest 
number of temporal derivatives and the windows with 
the steepest scalings with x, i.e. those that involve a; sin a; 
in our source convention. In principle we can also keep 
subleading terms in the analytic expression for the bis- 
pectrum. The rather involved results are neither illumi- 
nating nor necessary near the peak of the bispectrum as 
long as Xd 3> 1 and so we omit them here. 

It is straightforward to to derive an analytic expression 
for their contribution in the b — > limit. Integrating once 
by parts 



ho 

ho 
hz 



ds T)H - ci/2 



8 f 

ds o\ 



X(Ks), 



s f 



X(Ks), 



ds 
s 

ds 



c 2 s f 



f 



1 \ m + ai/2 

- v 1 ^— — X(Ks) 



X(Ks), 
(81) 



where X{x) — x 2 cosx. Combining these terms, we can 
approximate the bispectrum as 

G ^ A R (fc 1 )A R (fc 2 )A R (fc 3 ) 
hk 2 k 3 ~ 4^2 



(82) 



ds 
s 



clf + elf 



X(Ks) 



12 



The slow-roll parameters <j\ and t]h may be approxi- 
mated as [6] 



(l-c s )b 



dF 
din a' 



Cm 1 



2 1 



c Sl dF 
— b- 



(83) 



where recall that F is a step-like function in In a. As d — > 
0, dF/dhia can be approximated as a delta function with 
normalization set by J dlns(dF/d\na) ~ — f dF « 2. 
For finite step width d in field space, the inflaton traverses 
the step in As/s s w |dlns/d(/)|d. The window functions 
X oscillates on a time scale As = 1/k. Thus the integral 
is damped for ks s > <f>'/d and the same characteristic 
damping scale as in the power spectrum Xd = (f) 1 /nd (see 
Eq. 77) appears here. For the tanh step, the integral can 



be approximated following [5] 



kik 2 h 



A^(fci)A^(fc 2 )A^(fc 3 ) g (Ks\ X(Ks) 



2xdJ elf 



K 2 



2(:i -c,)~ -<■.. 



1 - Cs 
1 + c s 



where the damping function is 

T>{y) 



y 



sinh y 



(84) 



(85) 



In particular, for equilateral triangles, the bispectrum 
reaches a maximum of 



J.3 



10.77 



(2-c a )(l-c.) |6j 
12(1 + c.) C 2 



3(n s -l)/2 



at a scale of approximately fcs s w 22^- 

In Fig.[5]we compare this analytic expectation with the 
GSR and exact calculations at the peak of the equilateral 
bispectrum. For the analytic curves, we plot the analog of 



Eq. (84) for the terms in Eq. (81 ). For a small amplitude 



step, the agreement between all three is excellent. The 
i = 2 term (J20 + ^23) dominates at c s <C 1 and exhibits 
the cj 2 enhancement of the equilateral bispectrum over 
power spectrum effects which scale as 2(1 — c s )b/ (1 + c s ) 
through Eq. ( 76 1 6 . As c s — > 1 all terms are comparable 
and in fact partially cancel each other. While the total 
is suppressed in this limit, it scales in the same way with 
(1 — c s ) as the power spectrum. 

Fig. [3] verifies that the scaling remains qualitatively 
correct for the Planck amplitude step. The rough criteria 
for the validity of perturbation theory given in Eq. ( 79 ) 
becomes 

10-77 (2 ~ 2 (ff c ~ Ca) |^(2^) 3( "- 1)/2 A R < 1. (87) 

Thus sound speeds smaller than the c s = 0.15 example 
in the previous section would become allowed if y/bxd 
dropped significantly. 



We note that the result we find for the bispectrum here 
is significantly different from the estimate in Ref. |13) . 
The key difference is that Ref. [13] approximated the en- 
hancement to the bispectrum as arising for modes near 
horizon crossing due to the slow-roll parameter t]h be- 
coming large, and therefore the modes interacting more 
strongly near where they were freezing out. However, 
as we have demonstrated in this work, the key effect 
of the sharp feature in the background is that correla- 
tions are frozen at the time of the feature, s s , rather 
than near horizon crossing, s ~ 1/k. This means that 
correlations between curvature perturbations with mo- 
mentum k 3> l/s s are imprinted well before horizon 
crossing where the amplitude of the fluctuations is much 
larger and is oscillating with varying k. This results in a 
strongly scale dependent and oscillatory bispectrum that 
is only cut off at the damping scale derived above. In 
addition, the relative contribution of r]H and o\ terms as 
a function of width d was misestimated in Ref. [TBI. 



V. DISCUSSION 

In this work we have studied the bispectrum that arises 
in general models of single field inflation beyond the slow- 
roll expansion. Our integral approach allows the expan- 
sion history and inflaton sound speed to be arbitrary 
functions of time and encompasses all terms in the ef- 
fective field theory of inflation aside from those involving 
the extrinsic curvature and Galileon interactions. This 
form allows for a fast computation of all bispectrum con- 
figurations from a handful of one dimensional integrals 
and should facilitate efficient comparison of these models 
with data. 

The key assumption of our approach is that the expan- 
sion remains nearly de Sitter with ch <C 1 during infla- 
tion, while en is allowed to have arbitrary variation. In 
particular, we explicitly drop operators that contribute 
to the bispectrum as e 2 H and approximate the evaluation 
of the rest assuming that the curvature modefunctions 
remain perturbatively close to their de Sitter forms. 

Motivated by recent power spectrum analyses of 
WMAP [8j and Planck [S] data, we take as an illus- 
trative example a sharp step in the warp-brane tension 
in the context of DBI inflation jT3]. We show that our 
integral approximations are excellent for small amplitude 
steps in the warp, while for 40% steps remain accurate to 
20%. Further, we demonstrate that the consistency re- 
lation between the squeezed limit of the bispectrum and 
the power spectrum spectral index is respected in our 
approximations. 

A step in the sound speed due to a step in the warped- 
brane tension large and sharp enough to explain the 
high frequency oscillations in the Planck power spectrum 
would also generate a very large bispectrum peaking in 
the equilateral limit. For low sound speeds, c s , these bis- 
pectra are large enough to violate even weak criteria for 
the validity of perturbation theory and effectively put a 
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lower bound of c s > 0.15. Furthermore, as the sound 
speed increases, the step in the warp must become larger 
and larger in order to produce the same power spectrum 
features and ceases to remain viable as c s — > 1. In par- 
ticular, for c s > 0.5, one requires an order unity step in 
the warp to produce the favored 10% oscillations in the 
power spectrum. We have also demonstrated that other- 
wise degenerate scenarios that produce oscillations in the 
power spectrum - steps in the warp or potential - are dis- 
tinguished by their various bispectra. For c s < 0.5, the 
scenarios are clearly distinguished by the amplitude of 
their equilateral bispectra. 

Our technique also applies to the cases where the slow- 
roll approximation remains valid, but where slow-roll cor- 
rections may be large [ST]- We provide a complete and 
compact first-order expression for all bispectrum configu- 
rations. In Appendix |B 2| we demonstrate the consistency 
of our results with previous work |57j for several limiting 
cases. Given the observed value of the tilt, we show that 
large slow-roll corrections only exist in two cases: where 
one or more of the slow-roll parameters {e#, t]h , u\) is 
anomalously larger than the tilt n s — 1 or when the sound 
speed c s > 0.8. 

Our work is somewhat orthogonal to the work of Ref. 
|5"8] who evaluated the bispectrum for a class of P(X, <j>) 
theories allowing eu ~ 1 and <j\ ~ 1 while higher slow- 
roll parameters, <72, r/Hi etc., were assumed small. Fur- 
ther, Ref. |59j considered the bispectrum in Horndeski 
theories, relaxing the slow-roll assumption, again allow- 
ing for eji ~ 1 and <j\ ~ 1 in such a way as to preserve a 
scale invariant spectrum. However, in the region where 
our analyses overlap, one can demonstrate our results are 
equivalent to those of [55] (see Appendix B 2 and Refs. 

[SB ESI). 

The expressions presented here are ideally suited for 
use in concert with a fast estimator for the angular bis- 
pectrum in CMB. In particular this form should facilitate 
searches for the non-Gaussian counterparts to features in 
the power spectrum data which could confirm their pri- 
mordial origin. 
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Appendix A: Completing the Triangle 

In the main text, we deferred consideration of two op- 
erators that were not important in the DBI step example. 
In this Appendix we complete the GSR integral formu- 
lation for the bispectrum of general single field inflation. 



1. i = 3: 11* 

Beyond DBI models, general single field inflation mod- 
els can also have c~ 2 enhanced bispectra through the TZ 3 
term in the cubic action, Eq. (13). In the P(X, <j>) model 



it is given by Eq. (67). In the general effective theory of 



inflation it is associated with the M 3 mass scale defined 
in Ref. E21 or equivalently the C3 = — M|/M| coefficient 



of Ref. [60]. For a P(X, 

M A n {t)={-X) 



theory, in our notation 

Qnp 

X=X(t)' 



X(t) =- 



dt 



dX r 
2 



(Al) 



where <fi is the homogeneous background scalar field. 
With these conventions, M| = S(l — c^)/2 and 



Ml 
Ml 



1 A 



1-cjE 
1 



DBI 



(A2) 



(cf. [in] Eq. 5.12). For effective theories that parametrize 
turning trajectories in multifield models where the heavy 
degrees of freedom have been integrated out this scale 
becomes a specific function of the time-varying sound 
speed c s |19) . In summary, 





fo, 
1 


- 1 - 


[IJ 


(1 




(1 


2 C 2 

-k? 









t A 



= (!-<*)<*. 



DBI 

P(X, 

turn 
EFT 



(A3) 



It is interesting to note that 

T{cp) 
T{4>) - X 



P(X, 4>)=T{cp) 



1 - 



V&) (A4) 



would mimic the turning trajectory case at the effective 
field theory level. 

Inserting the 1Z 3 operator into the general expression 
Eq. j9j, we obtain 

n 1 !*!/" f ds 
— — / —S 3 (Ins) 



Q ^ J_ v , . 

hk 2 k 3 J AA 2 /CI ' 




(A5) 
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where 



aHsf 



(A6) 



As with the terms in the main text, an exact eval- 
uation of the 1Z 3 bispectrum contribution involves first 
solving numerically for the modefunction y through the 
exact equation of motion Eq. (27). Note that as x — > oo, 
y — > e ix and so the term in the second bracket tends 
to diverge as s 3 e~ 3tKs making it challenging to evaluate 
numerically. In practice when evaluating the exact bis- 
pectrum contributions, we regulate such expressions with 
an artificial damping factor at a sufficiently large s that 
the model is in the slow-roll regime. 

The GSR approximation can be constructed to avoid 
such problems. There are two types of terms in gen- 
eral: those that involve replacing y — ¥ yo and taking / w 
const, in Eq. (A5) and those that involve the first-order 



modefunction correction from Eq. (33) or /'//• The lat 



ter is required since S contributes at zeroth order in the 
slow-roll expansion. 
For the former case 



D 



A n {k 1 )An(k 2 )An(k3) 3 k x k 2 k 3 



A 2 



4 K 3 



ds 



S 3 (lns)(Ks) 3 sm(Ks), (A7) 



where we have replaced 1//* — > Ak as appropriate for ze- 
roth order modefunction expressions (see below). While 
compact in form, this expression is again numerically dif- 
ficult to evaluate at Ks — > oo. Instead, we first integrate 
this expression twice by parts and bring the result to the 



standard form of Eq. ( 36 ) 



uZul T&2 2^ T ^h 3 [.K), (A8) 



hk 2 k 3 " A A 2 . 

with the sources 



5*30 — 5*31 — S3, 
5*32 = S3, 



windows 



W 3 q(x) = x sinx, 

W 31 (x) = iy 32 (x) = cos x, (A10) 
and triangle weights 

T 3a =T 3 , T 31 = 3T 3 , T 32 = 2T 3 , (All) 

where 

3kik 2 k 3 



T 3 



K 3 



(A12) 



The modefunction correction terms involve both cor- 
rections to the external modefunctions evaluated at s* 



and nested integrals involving corrections inside the orig- 
inal integral in Eq. ( A5 ) . The former type is approxi- 
mated by the replacement of 1//* — » above. Note 
that the out of phase type contribution discussed in 
Ref. [5] vanish in this case since the (Ks) 3 cos(Ks) win- 
dow integrates to zero. The latter type involves integrals 
over (/(Ins) in Eq. (33). Part of the 1Z modefunction cor- 



rection term involves derivatives acting on / in Eq. ( A5 1 
This can also be brought into nested form by 



L 
f 



<iln s 



fJ 



(A13) 



The resulting terms can be simplified by integration by 
parts given that to leading order S 3 can be taken to be 
constant here. The result is that the new source becomes 



gS 3 — S 33 — S34 — S 35 , 



(A14) 



with now all terms combined into the form 

5 



G ^ A R (fc 1 )A R (fc 2 )A R (fc 3 ) 



4^ 

3=0 

+[T 3B I 35 (2k 3 ) + perm.]}. (A15) 
Here the additional windows are 
W33 (x) = x sin x + cos x, 
W 3i (x 



sin a; 
2 cos a;, 



, sin x cos x sin x 
W 35 (x) = 12( lcr -— 2 



4x 



(A16) 



and triangle weights 
3kik 2 k 3 



T 33 — 
T 34 = 
T35 = 



K 2 (K~2k 3 ) 
3hk 2 k 3 
~~ K 2 (K-2k 3 ) 2 
4fcifc 2 fc 3 



perm., 

[7(fci + k 2 ) — 3fc 3 ] + perm. 



(A9) T 3B 



K 2 (K-2k 3 ) 3 
-l-perm. 

2k]k 2 k 3 
~{K-2k 3 ) 3 ' 



[5(fci + k 2 f - 5(fci + k 2 )k 3 + 2k\] 



(A17) 



The final T 3 g term comes from the boundary term from 
unnesting the integrals through integration by parts. The 
boundary terms in from the modefunction terms can- 
cel since 



2 T 3i + (T 2 B+ perm.) = 0. 
3=3 



(A18) 



The Ssr boundary term in I 32 gives zeroth order slow- 
roll result for the P(X, <fi) model as 



3 k x k 2 k 3 ( 1 _ A 

„2 1 " 



hk 2 k 3 2 K 3 \c 2 ' "S 
which reproduces a well-known result [501 161] 



(A19) 
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2. i = 4: TLdOldiX 

Finally, for completenes s w e consider the (non-local) 
IZdilZdiX operator in Eq. ( 13 1 that was dropped in the 
main paper. Its contribution can be cast as a source 



Sa 



eg 

C 2 s f 



(A20) 



with windows 



Waq(x) = xsinx, 
Wa\{x) = cos IE, 



and triangle weights 



40 



4{K - 2k 1 ){K - 2k 2 ){K - 2k 3 ) 



kik 2 k 3 



(A22) 



T41 



1 



k 1 k 2 k 3 



K 



Note that in the squeezed limit of k\ = k$ <C &l — k 2 ~ 
&3, both weights scale as ks/kL and do not contribute. 
Hence they are not required for the establishment of the 
slow-roll consistency relation. On the other hand T41 
does contribute to the en terms in the slow-roll bispec- 
trum for other triangles. 



Appendix B: Slow-Roll Expansion 

With the completion of the GSR expression in <|A] we 
can expand the full result to first order in the slow-roll 
parameters e#, r\n and <j\ assuming they are nearly con- 
stant. This expansion trivially reproduces the usual slow- 
roll approximation for the i = 0,1,4 operators and im- 
proves their accuracy by defining the respective epochs of 
freezeout near horizon crossing. On the other hand, the 
i = 2,3 operators involve sources that are zeroth order 
in slow-roll parameters and thus complicate the slow-roll 
expansion [5U1 157] , 



1. GSR Derivation 

In the slow-roll limit, all of the sources in the GSR 
integrals can be taken to be slowly varying 



Sij(\ns) ps (In s* ) + 5- (Ins — Ins,), 



(Bl) 



where S'y is taken to be a constant. Thus the integrals 
reduce to 



hj = Sij (In s*)Wij(x*) + S-j 



dx 



Wij(x), (B2) 



where x* — Ks*. For the windows where linx^o Wtj 7^ 
0, the integral can be re-expressed as the evaluation of 
the source Sa 



SijQ.nsf)Wij(x*) 



(B3) 



WiAx) 



cos x e 
smx £l _ lE 



W(x/2) e 7/s -'< E 
W(x) e 7/3 - JE /2 
W n (x) e 2 -~< E /2 



TABLE IV. Freezeout epoch x ; for the various GSR windows 
where 7_b ~ 0.5772 is the Euler-Mascheroni constant, W is 



the power spectrum window of Eq. ( 38 1 and W n is the tilt 



(A21) window of Eq. (62 I 



at the freezeout epoch s/ = x//K where 

dx Wij(x) 



lncc/ — Ins, = / 



x Wij(x«)' 



(B4) 



In Tab. IV we give Xf for the various GSR windows. 



The GSR integrals can be further simplified by approx- 
imating the sources to first order in the slow-roll param- 
eters through 



aHs 



c s 
{aHs) 1 
aHs 

G' 



= — (lna)' ps 1 + e H + 01, 



2 /' 

-g w -2y ps 1 - n s 

(4e# -2t 1h + <j 1 ). 



(B5) 



With these approximations, the bispectrum to first order 
in the slow-roll parameters is given by 



1 - n. 



kik 2 k 3 



7ui + T\ )2 -!-:-! (^-1) T,„ -3ET.,, 



(Tn+T 12 ) 



aHs f J 4 



T 31 



tHrr (1 -c 2 s ){l~n s )+^ l rr 

4c 2 . 8 C 2 

Anik^Anik^Anih) 1 c s 



- 1 T 22 - ST f 



/ aHs 



t:i2 



(B6) 



where the final source is evaluated at Ks — e _7E and we 
have assumed the slow-roll parameters ch,t]h, and a% are 
constant in dropping terms. The accuracy can be further 
improved by evaluating the slow-roll parameters at their 
respective freezeout epochs rather than taking them to 
be constant. In fact the modefunction correction terms 
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contribute through these freezeout relations as 

7 

T2n = T23 + 2T25 + — T26 









T2B ( ln 2TT 


-0 


+ perm. 



^3n — T33 + 2T 34 + 2^35 









. 3B \ 2k 3 


-0 


+ perm. 



(B7) 



where we have used the identities Eq. ( 57 ) , ( A18 1 in elim- 
inating terms and the terms come from integrating 
the x sin x pieces of the windows by parts. Note that 
the power spectrum in the same approximation is 



G(ln*/) 



I+3K-I) 



(B8) 



where Sf = Xf/k — e 7 / 3 lE /2k. Eqn. (B6) reproduces 
and generalizes a well-known result for canonical (c s = 1, 
5 = 0) inflation g3 E] 



kik 2 k 3 fcifc 2 fc 3 



3 



(B9) 



Furthermore, for squeezed bispectra 

f 5 ^ g 

5 

= T-(4e H -2ry H + cr 1 )+C' 
12 \Kl 

as expected, while for equilateral bispectra 



108 cq _ 
35 /nl " 



24 

y 



A|(fc) 



(BIO) 



(Bll) 



1 



'A K (k)faHs 



1 - 



4 , 
21 ; 



2 / c s ST' 



. 7c 2 + 8 7c 2 - 10 



7 V«^s /7 ' ' 7c 2 " ' 7c 2 
2(c2-10)-135(c2-l)In(3/2) 



14c? 



(n s -l). 



The latter should be compared with the zeroth order 
in slow-roll result 



108 



1 



f cqO _ I 1 _ 

35 /nl ~ I c? ' 21" 



(B12) 



-l = -0.0674 




x Exact 
+ GSR 
— Analytic 
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FIG. 6. Slow roll corrections to the zeroth order equilateral 
bispectrum result f^. Agreement between the exact, GSR 
and analytic results is everywhere excellent for realistic values 
of the tilt n s — 1 = —0.0337. Shown also for comparison is a 
case of twice the observed value n s — 1 = —0.0674 showing that 
corrections are always of order the tilt. They only become 
fractionally large as c s — > 1 as the zeroth order result becomes 
suppressed. 



For example, for B = 0, we obtain using Eq. 

8c 2 - 17- 64(c 2 - l)ln(3/2) 



rcq 
J] 



NL 



pcqO 



pcqO 



7(cl 1) 
15e H + (Uje - 3)oi 
7(c* 1) 



(B13) 



Note that as one would expect, at small c 2 the correction 
is of order the deviation from scale invariance, n s — 1, 
unless en and G\ are anomalously larger. As c 2 — s- 1, 
the zeroth order term is itself suppressed and drops 
below the slow-roll corrections. Here the correction is 
large but the equilateral bispectrum itself is small. 



2. Examples 

We again illustrate our calculations with the DBI ex- 
ample, here with no feature (6 = 0). In Fig. [6j we show 
the exact, GSR and first order corrected analytic result of 



Eq. (B13) relative to the zeroth order of Eq. (B12l 



Compared with the zeroth order result, the corrections 
are of order the tilt n s — 1 throughout and are only frac- 
tionally large for c s > 0.8 where the zeroth order term 
itself drops below 0.2. The agreement between the exact, 
GSR and first order results are is excellent and consistent 
with deviations only at the next order in the slow-roll 
expansion (n s — l) 2 . By taking a second artificially high 
n s — 1, we also demonstrate that first order corrections 
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are always of order the tilt in these cases where the in- 
dividual slow-roll parameters are not anomalously larger 
than tilt. For both cases, most of the (n s — l) 2 deviation 
between the analytic and exact or GSR results can in fact 
be corrected by evaluating Eq. (Bll ) at the exact freeze- 



out epoch specified rather than using the approximate 
conversion in Eq. (B13). 



We can also compare the results derived here to known 
slow-roll results in the literature. In particular, we will 
compare to Ref. [57]. Since the expressions are rather 
large and cumbersome we will restrict again to the equi- 
lateral limit of DBI. Furthermore, the results of Ref. [57] 
are evaluated as an expansion about a fixed reference 
point. That is, while the dependence on wavenumber 
of coefficients in our Eq. (Bll) is implicit, their coeffi- 



cients are constant so that the resulting dependence on 
wavenumber is explicit. 

To compare results, we need to expand our expressions 
about some reference scale s*. Quantities evaluated at 
this time will be denoted with a subscript *. Working 
to linear order in slow-roll parameters and making use of 
Eq. JB5l, we find 



1 c s 



/ aHs 



1 



1 








C 2 " 




- n s * 


V 2 



l r 

7 



1 - — &U 
2(71* 



1 - c 2 
x ln(3fcs,e 7B ) 1, 



(B14) 



where Xf 
result, 

1 



Sks f = e lE . Finally, we make use of the 
1 2 - ln(2fcs*e 7B ; 



(1 - n s *) (B15) 



from Eq. ( |B8| ). Combining terms, for DBI inflation we 
find 



f<=q _ 
J NL — 



-{7(<4-l) + 15e ff * 



108c 2 , 

+ [8cl - 17 - 64(c 2 , - 1) ln(3/2)] (n sit - 1) 
+ [141n(3fcs,e 7B ) - 3]<t u |. (B16) 

In the limit where DBI inflation is described by power- 
law inflation at late times [10] . o\ = — 2e#, r\u = e#, 



and Eq. (|B16|) reduces to 
1 



JNL ~ 



35 
108 



1 - 



3-41n(3fcs*e 7B ) 



(B17) 



in agreement with Eq. (6.5) of Ref. [57 . We can also 
check the example of so-called generalized DBI inflation 
[52] . The slow-roll parameters there are defined as 



ev = 



1 (V_ 
V 



sgn^r 1 / 2 ) r 

W ~ V ' 3H T3/ 2 ' 



(B18) 



and are related to our slow-roll parameters by [57] 

e# = c s e v , rj H = - ~ey + ~r] V + jA 

3 

(Ti = -c s e v + c s n v - -A. (B19) 

Note that their r\ = 2(e# — T}h)- Working in the limit 
c s < 1, wc find 



/] 



35 



NL 



108c 2 , 



l + U^t [3 _14 In (3fcs,e 7E ) 



3A* 



+ [31 -I- 141n(3fcs*e 7B ) - 1281n(3/2)] 
2ey*c s * 



[43 - 7 In (3fcs,e 7B ) - 128 ln(3/2)] 
O(cl) (B20) 



in agreement with Eq. (6.11) of Ref. [57] . 



Appendix C: DBI Parameters 

In the main text, we choose the DBI parameters to 
have a fixed power spectrum amplitude A s and tilt n s — 1 
at s s for a model with no step 6 = 0. This condition fixes 
the parameters {</> e nd, As, /3} of the model and we set Vq 
to a constant as described in the main text. We relate the 
phenomenological and fundamental parameters through 
the slow-roll attractor 



H 



din a 



V2X 



V dlnV 



(CI) 



which determines X at <f>, and the definition of the sound 
speed c s = y/l — 2X/T. Through this relation, the phe- 
nomenological parameters p = {A Sl n s — l,c s } are de- 
fined given the field position which completes the set of 
fundamental parameters tz = {0, \g,/3}. Explicitly, we 
set 



,(tt) = y/1 - 2X/T 

V /dlnV 



-1/2 

+ 3T \ d<p J 
and use the Friedmann equation to set H 

3H 2 (tv) = (^-l) T{(j>)+V(ci>) 

The amplitude and tilt are then given as 

A s (n) 

(n s - l)(ir) 



(C2) 



(C3) 
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2nd(j)/dliia 
din A, dd> 



din a 



(C4) 
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This gives the phenomenological parameters as a nonlin- 
ear function of the fundamental parameters p(7r). To set 
the sound horizon at (f> to equal s s , we choose the appro- 
priate <^>end such that we effectively replace <f> with <^ e nd 
in the fundamental parameter set after the fact. 

Of course, we actually want the fundamental param- 
eters as a function of the phenomenological parameters 
7r(p). For c s <C 1, these relations are easily inverted 



For larger c s , the expressions are not readily invertible 
but from the solution at small c s , we can approximate 
small changes by linearizing and inverting the response. 
Starting from some parameter set po, we move to a new 
set p by iterating 



ln(7r/7r ) = 



<91np 
<9hi7r 



Mp/po), 



(C6) 




3?r V A s (n s - 1) ; 



(C5) 



until convergence. Here the inverse factor is the Jacobian 
matrix inverse. We repeat this procedure until we obtain 
all the desired values of the sound speed c s . 
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